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The relationship between pulsar-like compact stars and gravitational waves 
is briefly reviewed. Due to regular spins, pulsars could be useful tools for us 
to detect ~nano-Hz low-frequency gravitational waves by pulsar-timing array 
technique; besides, they would also be ~kilo-Hz high-frequency gravitational 
wave radiators because of their compactness. The wave strain of an isolate 
pulsar depends on the equation state of cold matter at supra-nuclear densities. 
Therefore, a real detection of gravitational wave should be very meaningful in 
gravity physics, micro-theory of elementary strong interaction, and astronomy. 
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1. Introduction 

Since the human mind first wakened from slumber, it has never ceased to 
feel the profound nature of space-time, especially the time-consciousness, 
in both philosophy and physics. However, a physical concrete of space- 
time is clarified only after Einstein's first insight: the space-time is a four 
dimensional continuum and the rule of the motion is the pure geometrical 
constrain that free particles follow the geodesies of the space-time, while 
the response of the space-time continuum to the matter is determined by 
the Einstein's equation in which a linear function of space-time curvature 
is in proportion to the energy-momentum of the matter. 

It is worth noting that the nature of space-time is a debating topic start- 
ing earlier than the general relativity and would not be terminated only by 
Einstein's pure geometrical arguments. Guided by different perceptions of 
space-time philosophies, many gravity theories are proposed, 1 with differ- 
ent interpretations of equivalence principle. Even in the general relativity, 
the equivalent principle is still a matter of debate. Therefore, experimental 
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tests of gravity theories, including in strong field and with fast motion, are 
critical to differentiate or falsify the gravity theories. Such experimental en- 
vironments are only available in astrophysics, especially related to compact 
stars known as white dwarfs, pulsars/neutron stars, and black holes. 

Topics of gravitational waves relevant to the pulsar astronomy are fo- 
cused on in this review. The binary pulsar tests for gravity theories are given 
in §2. Pulsars as tools of detecting and as sources of gravitational waves are 
presented in §3 and §4, respectively. Future prospects are discussed in §5. 

2. Pulsars, binary pulsars, and tests of gravity theories 

Soon after the discovery, pulsars are identified as a class of fast rotating 
rather than pulsating compact objects. It was used to believe that pulsars 
are neutron stars composed by hadronic matter, because of the very limited 
knowledge in sub-nucleon research at that time, but this view might not 
be true 2 since the lowest compact states could be of quark matters with 
strangeness rather than that of neutron liquid. Obscrvationally, there are 
two main categories of pulsars: the millisecond pulsars (MSPs) and the 
normal pulsars. Normal pulsars have rotational period from a few tens of 
milliseconds to a few seconds, while the MSPs' periods range from about 1 
millisecond to a few tens milliseconds. 3 Long term timing monitoring shows 
that the MSPs are more stable rotators compared to the normal pulsars, 
which might due to both observational reasons and pulsar intrinsic physics. 
For MSPs, the difference between model-expected time of arrival (TOA) 
of radio pulses and observed TOA is usually less than 10 percent of their 
periods, and the most stable MSPs can achieve ~ 100 ns level on the time 
scale of a few years. 4 

It turns out that most of the MSPs are in binary system. One particular 
interesting system is the recently discovered binary pulsar system, J0737- 
3039AB, where both of the two stars are radio pulsars. 5 The J0737 is also 
a highly relativistic celestial system. Binary pulsars with possible pulsar 
companions are listed in Table 1, obtained from ATNF pulsar catalog. 6 

Armed with such a kind of stable celestial clocks (i.e., pulsars) relativis- 
tically orbiting their companions, one can then test gravity theories in the 
case of strong gravitational fields, as illustrated in the classical system of 
PSR B1913+16. 7 In the J0737 system, two pulsars orbit each other with a 
period of 2.5 hours and a very low orbital eccentricity. Up to known, this 
J0737 system becomes the most relativistic binary system, the details of 
which can be found in the review by Kramer and Wex (2009). 8 

To test the gravity theories, one must compare the predicted TOA of 
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Tabic 1. Parameters for possible pulsar- neutron star systems 



Name 


P 


Pb 


a 


e 


J0737-3039A 


0.022699 


0.1023 


1.4150 


8.778c-02 


J0737-3039B 


2.773461 


0.1023 


1.5161 


8.778c-02 


J1518+4904 


0.040935 


8.6340 


20.0440 


2.495c-01 


B1534+12 


0.037904 


0.4207 


3.7295 


2.737c-01 


J1756-2251 


0.028462 


0.3196 


2.7564 


1.806c-01 


J1811-1736 


0.104182 


18.7792 


34.7827 


8.280c-01 


B1820-11 


0.279829 


357.7620 


200.6720 


7.946c-01 


J1829+2456 


0.041010 


1.1760 


7.2380 


1.391c-01 


J1906+0746 


0.144072 


0.1660 


1.4202 


8.530e-02 


B1913+16 


0.059030 


0.3230 


2.3418 


6.171c-01 


B2127+11C 


0.030529 


0.3353 


2.5185 


6.814c-01 



Note: P is the pulsar period in unit of second, Pf, is the orbit 
period in unit of days, a is the projected semi major axis in unit 
of light seconds, and e is the orbital eccentricity. 



a theoretical model with the observation. In this way, one needs thus the 
binary motion dynamical models, in which we put the gravity theory in. One 
can then calculate the theoretical pulse TOA at the solar barycenter. Note 
that, from a pulsar to the barycenter, various processes set in, including 
the photon propagation effects due to the gravitational field of both the 
binary system and solar system, the dispersion of pulsar radio signal due 
to interstellar medium and solar wind, and so on. One needs also the solar 
system ephemeris to convert the pulsar TOA at the radio telescope to the 
barycenter. The modeled TOA will then be compared with observed ones 
to see if the gravity theories is able to account for the observation. 

In reality, thanks to the phenomcnological framework of post-Kcplcrian 
(PK) parameters, with which gravity theories can be approximated, wc can 
independently measure these PK parameters by fitting the TOA data. A 
gravity test is then via checking the self-consistency of PK parameters for 
a particular gravity theory. There are 7 PK parameters which are possibly 
measurable in the near future: advance of periastron uj, gravitational red- 
shift parameter 7, Shapiro delay parameters r and s, orbit period derivative 
-Porb, spin-orbital coupling induced precession Q so . p and relativistic orbit de- 
formation 8g. These PK parameters, except 5g, have been measured in the 
0737 system. 

All the 7 parameters measured are functions of two unknown param- 
eters: pulsar masses, m a and m^. A double neutron star system is then 
overdetermincd if one detects three or more PK parameters. It is worth 
noting that the double pulsar system of J0737 offers an extra Keplcrian 
constrain, the mass ratio between two star, R(m a ,mb) = m a /mi, = Xb/x a , 
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with x the projected semimajor axes. Two recent reviews 8,9 are valuable in 
the topic of testing gravity with binary pulsars. 

3. Detecting gravitational waves with pulsars 

Directly detecting gravitational wave (GW) is the Holy Grail of present 
experimental researches, not only in gravity physics but also in astronomy. 
With the efforts since 1960s, 10 recent equipments (e.g., LIGO 11 ) may finally 
allow us to directly detect GWs although there is no confirmed detection 
now yet. In this section, we will review the ability of detecting gravitational 
waves using pulsar timing array (PTA). Potential roles of testing gravity 
with PTA are also presented here. 

GW is actually a perturbation of space-time, fully characterized by a 
wave-like metric perturbation. Detecting GW is thus identical to measure 
the wave-like metric perturbation 21 which can be performed by comparing 
geodesies of two test objects approaching to and departing from each other. 
Such experiments fall into four categories: 1. Tracing the motion of two free- 
falling test objects (e.g. LIGO, LISA, GEO, TAMA, and so on), 2. Detecting 
the deformation of finite extend solid body (e.g. Bar detector, Sphere detec- 
tor, and so on), 3. Measuring the Doppler shift of electromagnetic signals 
from distance free-falling objects (e.g. Doppler tracking of satellite, pulsar 
timing array, laser ranging, LISA), 4. Checking the perturbation of a cos- 
mological system (e.g. cosmic background B mode detection, weak lensing 
survey). Among all these possible ways, PTA is one of the promising tech- 
niques to directly detect gravitational waves, being unique to detect GW 
at nano-Hertz band. 12 

As we have shown, MSPs are very stable celestial clock in the Galaxy. 
GWs perturb the background space-time of the Galaxy, such that pulsar 
pulse signals get red or blue Doppler shift along the path from pulsar to 
earth. It turns out that such GW-induced frequency shift only involves the 
metric perturbation at the pulsar and that at the earth. The GW-induced 
frequency shift can be obtained to be 13 



a It should be born in mind that detecting of GW is not detecting any types of metric 
perturbation. GW detection focuses on detecting the oscillatory part of the metric per- 
turbation with strain h decrease as r _1 , such that gravitational wave could carry energy 
and momentum to the infinity. 



[hijit^-hijit-D/cB)}, 



(1) 



to 



2 (1 + n g • n) 
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where w is the pulsar angular frequency of spin, n is the pulsar direction, 
D = Dn with D the distance to the pulsar, n g is the GW propagation 
direction, and h is the perturbation of metric. The GW-induced timing 
residuals in pulsar TOA is therefore R = J Au/ujdt. 

Due to the intrinsic noises and possible non-modeled accelerations of 
pulsars, it is unlikely that one can use R(t) of a single pulsar to detect GWs. 
Nevertheless, magic happens if we correlate the residuals of Rj and Rj of 
two pulsars. From Eq. (1), one can have a correlation of (RiRj) = C(9)a 2 
for two different pulsars in general relativity, with a the RMS (root-mean- 
square) of a single pulsar's residual. Note that the correlation C(9) is a 
determined function only involving the angular, 9, between two pulsars, 14 
and this correlation C (9) certainly plays a vital role in detecting GW using 
a array of pulsar timing data (PTA) since the shape of C{9) is uniquely 
determined by a gravity theory and there is no other physical processes 
to make the pulsar signal correlated for two pulsars widely separated with 
a distance of several thousand light years away from each othcr b . In the 
general relativity theory of gravity, fortunately, the correlation C{9) has a 
very simple form of 13,14 

CW = ^-f + i(l + ^)), (2) 

where x = (1 — cos#)/2. 

We may make sense of the C(#)-curves from simple symmetric reasons. 
If a monochromatic general relativistic GW is propagandizing along 'z- 
axis' direction (there will be 180° symmetry and 90° anti-symmetry in x-y 
plane), then correlation C{9) between two pulsars with 0° or 180° angular 
separation are positive, while C{9) for 90° will be negative. This make a 
U-shaped C(6>). 14 Note that C(180°) ^ C(0) which will be explained later. 
One can then measure such multi-pulsar correlation to detect GWs. Jcnct 
et al. (2005) 12 had investigated the statistical properties of such detection 
processes. Their results show that regular timing observations of 40 pulsars 
each with a timing accuracy of 100 ns will be able to make a direct detection 
of the predicted stochastic background from coalescing black holes within 
5 years. We compare the detection abilities for GW detectors in Fig. 1 , for 
GW background due to coalescing supermassive binary black holes (BBH). 

From symmetry arguments above, it is clear that the shape of C{9) de- 



b The imperfectness of terrestrial clock and un-modeled solar system dynamics may intro- 
duce also correlation between measured Ri of pulsars, however the angular dependence 
of such correlations is very different from that C(6) presented in Eq.(2). 
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Fig. 1. A comparative for abilities of detecting GWs for various GW detectors, where 
the x-axis and y-axis are the GW frequency and the characteristic strain of GW respec- 
tively. Most of the labeling are self-explanatory. The PTA curves are at in the left-top 
corner, where the length of data and pulsar average timing noise level are also labeled. 



pends on the polarization of GW (see Fig. 2). Einsteins theory of gravity 
predicts waves of the distortion of space-time with two degrees of polariza- 
tion; alternative theories predict more polarizations, up to a maximum of 
six. 15 Lee et al. (2008) 13 analyzed such polarization effects and conclude 
that for biweekly observations made for five years with rms timing accu- 
racy of 100 ns, detecting non-Einsteinian modes will require: 60 pulsars in 
the case of the longitudinal mode; 60 for the two spin-1 'shear' modes; and 
40 for the spin-0 'breathing' mode. Further more, they showed that one 
can test gravity theories by checking GW polarization, i.e., to discriminate 
non-Einsteinian modes from Einsteinian modes, we need 40 pulsars for the 
breathing mode, 100 for the longitudinal mode, and 500 for the shear mode. 
These requirement is beyond present observation technology, but could be 
easily achieved using SKA or FAST telescope. 16 ' 17 

Another interesting topic on detecting GWs using PTA is about the 
dispersion relation of GW, 18 since the function of C{9) and the detection 
statistics depends also the mass of graviton. It is found 18 that C(180°) 
increases to match the value of C(0) as the graviton mass increases (see 
Fig. 3 for details). In the case of massless GW background, we know that 
the GW has 180° degree symmetry due to the polarization property, but 
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Fig. 2. The C{9) curves for different kinds of GW polarization, with a denoting the 
power index of the GW background. 



why C(0°) ^ C(180°)? It turns out that GW propagation breaks up this 
180° symmetry by the geometric factor in Eq. (1), which reads 1 + e z ■ n for 
the massless case. For the case of a massive GW background, the geometric 
factor reads 1 + — n K • n, where the graviton mass reduces the asymmetry. 
For the limiting case, where the GW frequency is just at the cut-off fre- 
quency, the dispersion relation tells us that such a GW is not propagative, 
then the 180° symmetry is restored. Therefore, we would expect that the 
correlation function are of 180° symmetry for very massive gravitons. 

Lee et al. (2009) further find that it is possible to measure graviton mass 
using PTA and one will get 90% probability to differentiate between the 
results for massless graviton and that for graviton heavier than 3 x 10~ 22 
cV, if biweekly observation of 60 pulsars are performed for 5 years with 
pulsar RMS timing accuracy of 100 ns in the future. 

As we have shown that PTA can be constructed to measure the alter- 
native polarization modes of GW and the GW dispersion relation. These 
measurements provide tests for gravity theory in the weak field/high veloc- 
ity region, which are different from that of the solar system tests (i.e., the 
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weak field/slow velocity case) and the binary pulsar tests (i.e., the strong 
field/slow velocity cases), because it is not completed to describing GW us- 
ing post-Newtonian formalism and the scalar and tensor sectors of gravity 
theories arc different. 19 



4. Gravitational wave radiation from pulsars 

Pulsars are not only as tools to detect GWs, but also strong GW sources 
because of their compactness and the rapid mass changes. Indirect evidence 
for GW from binary pulsars has been discussed in §2, whereas a direct detec- 
tion of GW from pulsars with ground-based facilities should be meaningful. 
It is recognized that the GW amplitudes of isolate pulsars depend on the 
equation of state (EoS) of cold matter at supra-nuclear density, which is 
strongly related to the understanding of QCD (quantum chromo-dynamics) 
at low energy scale, still another challenge for physicists today. A mixture 
of quantum (QCD) and gravity (relativity) makes this project more funny. 

We are still not sure about the nature of pulsar-like compact stars 
though discovered since 1967. It is conventionally believed that these com- 
pact stars are normal neutron stars composed of hadronic matter, but one 
can not rule out the possibility that they are actually quark stars of quark 
matter 2 (see, e.g., a review 20 ). Quark stars with strangeness are popu- 
larly discussed in literatures, which are called as strange (quark) stars. The 
EoS of realistic quark matter in compact stars, depending non-perturbative 
QCD, was supposed to be of Fermi gas or liquid, but could be of classical 
solid in order to understand different manifestations of pulsar-like stars. 21,22 

Besides QCD, that pulsars are quark stars should also be meaningful in 
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GW physics. 23 (i). GW being EoS- dependent. Rotation (r) mode instability, 
which would result in GW radiation, may occur in fluid quark stars if the 
bulk and shear viscosities of quark matter is not sufficiently high, but no 
r-mode instability occurs in solid quark stars. Even in case of solid quark 
stars, the GW amplitude is relevant to the quadrupole deformations 27 (e.g., 
mountain building on stellar surface) sustained by elastic or magnetic forces 
on stellar surface. A quark or neutron star with quadrupole deformation 
would be a GW radiator if it has precession either free or torqued, and the 
precession amplitude (or the angle between spin axis and spindle of inertia 
ellipsoid) is determined by EoS 24 and determines GW strain, (ii). GW being 
mass-dependent. A very difference between quark and normal neutron stars 
is that the latter is gravity-bound while the former is confined additionally 
by self strong interaction, that results in the fact that quark stars could be 
very low massive 25 (even to be of ~ 1O _3 M ) but neutron star cannot. Low 
mass quark stars, either in liquid or solid states, are surely very weak GW 
emitters. This mass-depend nature makes it more complex to constrain EoS 
of target compact stars by negative results of LIGO GW detections. The 
points of above (i) and (ii) are certainly very useful for us to observationally 
distinguish quark stars from normal neutron stars in the future. 

Pulsars spin usually at frequencies > 10° Hz, and we thus are interested 
in LIGO to detect their GWs, from Fig. 1. There are two kinds of GWs 
from pulsar-like compact stars: continuous GWs due to spin and bursting 
GWs due to stellar catastrophic events (e.g., star quake 26 or binary coales- 
cence). It is worth noting that all the upper limits estimated from LIGO 
science runs depend on simulated waveform types of GWs (i.e., astrophysi- 
cal GW radiative mechanisms). For continuous GWs, the waveforms could 
be better understood, and their searches are significantly more sensitive, 
especially when informed by observational photon astronomy and theoret- 
ical astrophysics. 28 The waveform of bursting GW is a matter of debate, 29 
and such kind of GW searches is also focused by LIGO, especially on the 
super-flares of soft gamma-ray repeaters. 30 

5. Summary and Future prospects 

Pulsars could be useful tools to detect GWs by PTA technique, they would 
also be strong GW radiators; a real detection of gravitational wave should 
be very meaningful in gravity physics, micro-theory of elementary strong 
interaction, and astronomy. A successful detection of GW by PTA may 
provide a test of GW polarization and measure the graviton mass. Thought 
indirect evidence for GWs from pulsar timing in binaries has been obtained, 
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a direct GW detection of pulsar-like stars is also expected as persistent or 
transient sources. The strain of GW from an isolated compact star depends 
on the equation of state of cold matter at supra-nuclear densities. 

Pulsar timing array projects are promising for detect GWs. We can 
achieve h c = 10~ 15 region for several year continuous pulsar timing moni- 
toring. If bi-weekly observations are made for five years with RMS timing 
accuracy of 100 ns, then 40 pulsars are required for general relativistic 
modes, 60 for the longitudinal mode; 60 for the two spin-1 shear modes; 
and 40 for the spin breathing mode. Additionally, we may measure the 
graviton mass through PTA techniques. With a 5-year observation of 100 or 
300 pulsars, we can detect the graviton mass being higher than 2.5 x 10~ 22 
and 10 -22 eV, respectively. Ultimately, a 10-year observation of 300 pulsars 
allows us to probe the graviton mass at a level of 3 x 10~ 23 eV. 

For the task of measuring the GW polarization and the graviton mass, 
there is one critical requirement: a large sample of stable pulsars. Thus the 
on-going and coming projects like the Parkes PTA, 31 the European PTA, 32 
the Large European Array for Pulsars, 33 the FAST 16,17 and the SKA would 
offer unique opportunities to detect the GW background and to probe into 
the nature of GWs, both physical (the GW polarization and the graviton 
mass) and astronomy (the GW sources). 

Other important requirements for a successful PTA include high stabil- 
ity of pulsar intrinsic noises and a low measurement noise. We need pulsar 
survey with better sky coverage as well as good observing system, especially 
with better band width 34 to get better signal to noise ratio and to subtract 
the interstellar medium effects. Better radio frequency interference filtering 
technology will also be very helpful such that we can use the full band data 
and reduce the terrestrial contamination. Better timing techniques (such as 
timing in full Stokes parameters 35 ) could also be preferred. 

Acknowledgments 

We would like to acknowledge useful discussions at the pulsar group of 
PKU. This work is supported by NSFC (10833003, 10973002), the Na- 
tional Basic Research Program of China (grant 2009CB824800) and LCWR 
(LHXZ200602). 

References 

1. C. Will, Living Reviews in Relativity 9, 3 (2006). 

2. J. M. Lattimer, M. Prakash, Science 304, 536 (2004). 

3. D. R. Lorimer, Living Reviews in Relativity 8, 7 (2005). 



September 14, 2011 0:15 



WSPC - Proceedings Trim Size: 9in x 6in ms 



11 

4. J. P. W. Verbiest, M. Bailes, W. A. Coles, et. al., MNRAS 400, 951 (2009). 

5. A. G. Lyne, M. Burgay, M. Kramer, et al., Science 303, 1153 (2004). 

6. R. N. Manchester, G. B. Hobbs, A. Teoh, et al., AJ 129, 1993 (2005). 

7. T. Damour and J. H. Taylor, Phy. Rev. D 45, 1840 (1992). 

8. M. Kramer and N. Wex, Classical and Quantum Gravity 26, 073001 (2009). 

9. I. H. Stairs, Living Reviews in Relativity 6, 5 (2003). 

10. R. L. Forward, D. Zipoy and J. Weber, Nature 189, 473 (1961). 

11. B. P. Abbott, et al., Reports on Progress in Physics 72, 076901 (2009). 

12. F. A. Jenet, G. B. Hobbs, K. J. Lee, et al., Ap J 625, L123 (2005). 

13. K. J. Lee, F. A. Jenet and R. H. Price, ApJ 685, 1304 (2008). 

14. R. W. Hellings and G. S. Downs, ApJ 265, L39 (1983). 

15. D. M. Eardley, D. L. Lee and A. P. Lightman, Phys. Rev. D 8, 3308 (1973). 

16. R. Nan, Q. Wang, L. Zhu, et al., CJAA S6, 020000 (2006). 

17. R. Smits, D. R. Lorimer, M. Kramer, et al., A&A 505, 919 (2009). 

18. K. J. Lee, et al. (2010), in prepration. 

19. M. Maggiore, Gravitational waves (Oxford: Oxford University Press, 2008). 

20. R. X. Xu, J. Phys. G: Nucl. Part. Phys. 36, 064010 (2009). 

21. R. X. Xu, ApJ, 596 L59 (2003). 

22. R. X. Xu, in Compact stars in the QCD phase diagram II, May 20-24, 2009, 
KIAA-PKU, Beijing (arXiv:0912.0349). 

23. R. X. Xu, Astroparticle Physics, 25, 212 (2006). 

24. I. H. Stairs, A. G. Lyne, S. L. Shemar, Nature, 406 484 (2000). 

25. R. X. Xu, Mon. Not. Roy. Astron. Soc, 356, 359 (2005). 

26. R. X. Xu, D. J. Tao, Y. Yang, Mon. Not. Roy. Astron. Soc. 373, L85 (2006). 

27. B. J. Owen, Phys. Rev. Lett., 95 211101 (2005). 

28. B. J. Owen, preprint (arXiv:0904.4848). 

29. P. Kalmus, PhD thesis submitted to Columbia Univeristy (arXiv:0904.4848). 

30. J. Horvath, Mod. Phys. Lett. A20 2799 (2005). 

31. G. B. Hobbs, M. Bailes, N. D. R. Bhat, et al.„ PAS A 26, 103 (2009). 

32. B. W. Stappers, M. Kramer, A. G. Lyne, et al., CJAA S6, 020000 (2006). 

33. B. Stappers, W. Vlcmmings and M. Kramer, in Proceedings of the 8th Inter- 
national e-VLBI Workshop. 22-26 June 2009. Madrid, Spam (2009). 

34. X. P. You, G. Hobbs, W. A. Coles, et al., MNRAS 378, 493 (2007). 

35. W. van Straten, ApJ 642, 1004 (2006). 



